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On-axis gain of antennas with rough reflecting surfaces is computed 
as a junction of rms surface deviation e, correlation distance c, antenna 
area A and wavelength X. Gaussian stationary surface deviations, Gaussian 
correlation functions, and uniform illumination are assumed. Antennas 
with rectangular and circular apertures are considered. It is shown that 
a normalized gain can be defined which has the same functional form for 
both. A principal result of this work is quantitative calculation of the 
on-axis antenna gain when the normalized variance (Awe/X) of the rough 
surface is larger than 4. The off-axis gain is also considered, and it is 
shown that in the asymptotic limit (as X — » 0), the gain reduces to that 
obtained by using geometrical optics. 

I. INTRODUCTION 

The gain of shallow paraboloid reflector antennas with random 
surface deviations has been derived by Ruze. 1 ' 2 The deviation was based 
on scalar Kirchhoff approximation to the radiation from reflector 
antennas. The surface deviations were assumed to be gaussian stationary 
with gaussian correlation functions. On these bases, an approximate 
solution for the antenna gain was obtained in terms of an infinite 
series. The series has been evaluated for relatively small rms surface 
deviations, e, in comparison to the wavelength, X, namely (47re/X) ^ 4. 
Asymptotic limits (as X — > 0) for the gain were also given by Ruze 2 
based on a similar analysis by Scheffler 3 . On-axis gain measurements 
of large reflector antennas as a function of frequency, exhibit the char- 
acteristics as predicted theoretically by Ruze. 

The present work was motivated primarily to determine the gain 
in the intermediate region between very long and very short wave- 
lengths and to establish a criterion for applicability of the asymptotic 
limit. Of primary interest was the near axis field distribution in the 
focal plane of a paraboloid reflector antenna illuminated by an inci- 
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dent plane wave. However, since both the far-field radiation pattern 
and the field distribution in the focal plane are Fourier transforms 
of the antenna aperture illumination, the derivations by Ruze are 
applicable for determining both the far-field and focal-plane dis- 
tributions. 

The series solution for the antenna gain obtained by Ruze does 
not seem to be suitable for numerical computations for large argu- 
ments. This is because some of the terms in the series will assume 
large values before the terms begin to decrease. However, the series 
for the on-axis gain is related to an exponential integral. The expo- 
nential integral also has an asymptotic series representation, which 
is particularly suitable for numerical computation for large argu- 
ments. On this basis, the on-axis gain has been computed as a func- 
tion of the rms surface deviation to the wavelength ratio and for a 
range of correlation parameters. The asymptotic limit for the gain 
is evident from these computations. 

The off-axis gain is also considered. Asymptotic representations of 
the series which may facilitate the off-axis gain computations are dis- 
cussed. The limiting value (A -* 0) for the off-axis gain is obtained, 
and it is shown that in this limit, the gain reduces to that obtained 
from geometrical optics.* 

The gain of antennas with rectangular apertures and gaussian 
stationary surface deviations is presented by assuming uniform il- 
lumination. A generalization to include certain types of nonuniform 
illuminations is discussed. The on-axis gain for antennas with cir- 
cular apertures also is given. It is shown that the on-axis gain for 
antennas with rectangular and circular apertures can be normalized, 
such that the normalized gain is the same for both. The off-axis gain 
is expressed in terms of series with known asymptotic expansions. 

II. ANTENNA GAIN 

The far field gain, G(0,$), in the vicinity of the axis of a shallow 
paraboloid reflector antenna with surface deviations, z(x,y) is, using 
the scalar Kirchhoff approximation 5 

G(0, ft) = ^ 

fj ff E„(x, y)E*{x l , Vl ) exp (j{0 x u+Pj)+2k[z(x, y)-z(x i , ?/,)]!) dsds, 

fj E a (x, y)E*(x, y) ds (1) 
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where E a is the projected electric field on the antenna aperture and s 
is the aperture area. 

j. _ £* _ free space propogation constant 
A 

A = wavelength 

p t = k sin cos $ (2) 

ft, = /c sin 6 sin $ (3) 

and $ are the spherical coordinates indicated in Fig. 1. 

u = x — Xi (4) 

v = y - Vi (5) 

The Kirchhoff approximation is based on the assumption that the 
surface is locally plane, and hence equation (1) is applicable to sur- 
faces for which the curvatures are small. 

Equation (1) can also be used to determine the power distribution 
in the focal plane of shallow paraboloid reflector antennas in the 
vicinity of the focal point, in which case (referring to Fig. 1) 

0. = fee/// (°) 

A, - ky,/f (7) 

where x, and y, arc the coordinates in the focal plane and / is the 
focal length. 

If z(x,y) is a Gaussian stationary random variable with zero mean 
it has been shown 1 - c that by performing the statistical averaging, the 
expectation value for the gain, {G{0, *)) llv is: 

(G(6, $)) av = ^exp(-5 2 ) 

[[ [[ E a (x, y)E*{x i , yd exp Wm + M] exp [5 2 r(w, »)] rfs tfs, 
•/J. JJ. (8) 

E a (x, y)E* a {x, y) ds 



// 



where 



8 = ^e (9) 

A 
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z(x,y) 




Fig. 1 — Antenna coordinates. 

e = rms surface deviation 
b'r{u, v) = correlation function. 

To evaluate equation (8), four integrations have to be performed. 
It is shown in Appendix A that for antennas with rectangular aper- 
tures two integrations can be readily eliminated for certain types 
of illuminations, truncated cosine illuminations, for example. 

In particular for uniform illuminations, E„(x, y) = 1, and for a 
gaussian correlation function with 



r(u, v) = exp ( — 



w" + v'- 



(10) 



where c is the correlation length, it is shown in Appendix A that the 
expectation value of the gain is: 



<0(0,»)>«r = CXp(-5 2 X? o (0, *) 



- A„ 



(11) 



where Go(d, $) is the antenna gain in the absence of surface devia- 
tions. 
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For an antenna with a rectangular aperture 

n (a *> 4tt,-1 / sin fl y a sin v b \ 2 n ~ 

where ^4 = 4ab is the aperture area. 



= ^sm0 (13) 

A 



and 

A„ < 



K«n) f 



H. 



(14) 



Equation (11) agrees with the gain derived by Ruze for antennas 
with circular apertures except for the term A„. This term is small, if 
the correlation distance c is small compared with the linear dimen- 
sions of the antenna. This assumption is made in the subsequent 
computations. 

For antennas with circular apertures the exact evaluation of equa- 
tion (8) is in general more difficult. However, for uniform illumination, 
the on-axis gain ((7(0, 0))„ v is evaluated exactly in Appendix B with 
the aid of Q functions. The gain has the same functional form as equa- 
tion (11) with p = 0. In particular for n/2(D/c) 2 » 1 

2l ' (15) 



" DdrnY 
where D is the antenna aperture diameter. 

III. ON-AXIS GAIN 

Equation (11) can be readily computed for small values of 8 2 . For 
large values of S 2 the terms h 2n /n\n will become very large; therefore, 
the series is not suitable for direct computation if 8- is large. Never- 
theless, the gain on-axis can be readily computed by noticing that 
the scries in equation (11) for f3 = is related to an exponential 
integral, which also has an asymptotic representation. 

The exponential integral, E t , can be written 7 

E<(x) -7+Inx+E-^; (16) 

where y is Euler's constant. The asymptotic series {x -> oo) for 
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Et(x) is 

Though the asymptotic series diverges for all finite values of x it can 
be used to evaluate Ei(x) for large x by using up to N terms, 8 where 
N is an integer nearest to the value of x. 

In terms of the exponential integral, the on-axis gain for both 
rectangular and circular aperture antennas is 

«?«), 0)),, = (&)' 

•js 2 exp (-«*) + (~-J S 2 cxp (-r')[ff,(5 2 ) - In 5 8 - 7]} (18) 
where D is related to the antenna area, A, by 

4 = £|5. (19) 

One parameter in (18) is readily eliminated by defining a normal- 
ized on-axis gain, (g (0,0) ) nv , by 



(a(Q 0)) - @& o»~ 



= 8 2 exp (-5 a )|l + [§J EBiW - ln *" - T]}- 



(20) 



The normalized gain thus depends only on two parameters, 5 2 and 

(c/D y. _ 

Equation (20) has been computed by using a share program for 
the computation of the exponential integral*. This program computes 
E { (x) with at least four-decimal accuracy. 

Computations have been performed for 10~ 4 =g 5 2 ^ 80 and for 
10~ 3 ^ c/D ^0.1. The computed normalized gain is shown in Fig. 2. 

The computations show the normalized antenna gain has three 
distinct regions which are characterized by the normalized rms surface 
deviation to wavelength ratio, 5. 

In the region ^ 8 2 ^ 1 the normalized antenna gain is nearly 
independent of the correlation length c, and increases almost linearly 
with 5 2 . In the region 1 ^ 5 2 ^ 20 the gain is dependent on both S 2 

* Contributed by D. S. Villars. 
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Fig. 2 — Normalized antenna gain. 

and c. In the region 8~ > 20 the gain is almost independent of 5 2 and 
is a function of c/D only. This region is the asymptotic region. For 
the range of parameters used in the computation, the gain in the 
asymptotic region, for a given c/D , deviates by less than 5 percent 
from the asymptotic value. 

The curves shown in Fig. 2 seem to confirm the general charac- 
teristics of the measured gain as a function of frequency of large 
reflector antennas presented by Ruze. 2 The presented measurements 
extend only slightly from the first into the second region but not 
sufficiently far to determine qualitative agreement between the theory 
and experiment in much of the intermediate and all of the asymptotic 
regions. A detailed comparison of the measured and computed gain 
can not be made since uniform illumination has been assumed in the 
computation. 

IV. ASYMPTOTIC VALUES FOR THE OFF-AXIS GAIN 

Computation of the off- axis gain directly from equation (11) can 
only be readily performed for relatively small values of 8 2 . An 
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alternate representation for the off-axis gain for large values of 8 
is obtained by expanding the exponential in the second term of equa- 
tion (11) in a power series. By using this expansion and neglecting 
A„, (11) can be rewritten as follows: 

(G(e, ft)) av = exp (-8 9 )G (e, ft) + (|) 2 ± (-1)'" ^p (I) 2 '" (21) 

where 

oo / «2\ n + m + 1 

A m (8 2 ) = ex P (-5 2 ) S ;!, n -+» • (22) 

The series for A m {l-) are special cases of functions considered by 
Barnes 9 who obtained their asymptotic expansions. These functions 
were also studied by Ford 10 who also presented a recurrence relation 
for the coefficients of the asymptotic series. The above functions 
designated by Gp(x,®) are: 

<?,(*, ®)= t - *" • (23) 

n=o n\ (n + ©) 

The functions A OT (8 2 ) can be written as: 

-*«(*■) - (5 2 ) m+2 exp (-fi a )(? M+a (5 a , 1). (24) 

Only the asymptotic limit for the off-axis gain is considered. For 

x -» oo 



(?,(*, 9) = S*& [l + 0@] 



hence for 8- -> co 

^, n = l + o(^- (25) 

The off-axis asymptotic gain will be designated by G{6, <J>) M and is 
given by: 

G(9, •). = (^) 2 exp [ - (^ sin e)j (26) 

The corresponding normalized gain is found as 

g(e, •). = (I;) 2 exp -(gam «) ]■ (27) 



ANTENNA SURFACE DEVIATIONS 1645 

The asymptotic value for the gain, equation (26), is in agreement 
with the gain obtained, based on an approximation to the gaussian 
correlation functions, 3 and also with results obtained by using geome- 
trical optics. 4 Equation (24) is independent of frequency but strongly 
dependent on the ratio e/c. This ratio has been interpreted as an 
average surface slope. 2 The range of 8 over which the above off-axis 
approximation applies has not been determined precisely, however, 
it is reasonable to assume that this range will correspond to the asymp- 
totic region for the on-axis gain. 

V. CONCLUSIONS 

The on-axis gain of antennas with gaussian stationary random 
surface deviations and gaussian correlation functions has been deter- 
mined for antennas with rectangular and circular apertures by assum- 
ing uniform illumination. For both types a normalized expression for 
the gain was derived which depends only on the normalized rms sur- 
face deviation e, to wavelength A ratio, 8(=47te/a), and the ratio of 
the correlation length c to a defined linear antenna dimension D . 
For circular antennas, D is the diameter. 

The antenna gain as a function of 5 exhibits three distinct regions: 
(*") ^ S 2 ^ 1, (it) 1 ^ 3 2 ^ 20, and (m) S 2 > 20. The last is called 
the asymptotic region. In this region the gain is nearly independent of 
wavelength. 

The computed gain exhibits in general the characteristics of the 
measured gain as a function of frequency of large reflector antennas 
reported in the literature. These measurements extend only partially 
into the second region and have not been obtained in the third 
(asymptotic) region. 

For large values of S 2 the off-axis gain can be expressed in terms 
of series with known asymptotic expansions. The limiting value for 
the off-axis gain has been obtained and reduces to that obtained by 
geometrical optics. 

VI. ACKNOWLEDGMENTS 

The author wishes to thank B. H. Bharucha for the valuable dis- 
cussions on the properties of the Q functions, and H. G. Cooper for 
helpful comments. J. A. Arnaud brought Ref. 4 to the author's 
attention. 



1646 



THE BELL SYSTEM TECHNICAL JOURNAL, OCTOBER 1968 



APPENDIX A 



Gain of Antennas with Rectangular Apertures 

To evaluate equation (8) for antennas with rectangular apertures, 
consider the following integral 

</> ftV - exp (-« a ) J [J f E &* y) B <&* . yJ 

•exp [8 2 r(u, v)] exp \j(fi,u + /3„v)] dx dxj dy dy x (28) 



with 



u = x — x. 



v = y - yi. 



(29) 
(30) 



Since equation (28) contains the correlation function in terms of 
u and v, it is preferable to introduce the u,v coordinate system. 

In the x,x x coordinate system the integrations are over the square 
region shown in Fig. 3a. In the y, y x system the region is similar. With 



(-a,a) 




(-a,-a) 



(a,-a) 



(a,-a) 



fa) 




(~a,a) 



Fig. 3 — Coordinate transformation. 

the coordinate transformation equation (29), the transformed region 
in the x lt u coordinate system is also shown in Fig. 3b. 

In the x lf u plane the integration with respect to .r x is readily per- 
formed for certain types of illumination functions.* In particular, let 



E a (x, y) = E ax (x)E au (y) 



(31) 



* A similar method has been used by Hoffman in his treatment of scattering 
of electromagnetic waves from a random surface. 6 
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where 

E„(xi + u) = E g n (xt)U(u) (32) 

n-l 

with a similar equation for E ay {yi + v). An example of such illumina- 
tions are truncated cosine illuminations, where equation (32) will 
consist of two terms. 

It is sufficient to consider the following integral 

</i>- = /" /" ' g(x l )Ku,v)dudx i . (33) 

Referring to Fig. 3, equation (33) can be written as 

COuv = J J gfaMu, v) dXi du 

+ f f gixJKu, v) dx, du (34) 
let 

0fe) = / g(xO dx } 
then 

</,>„ = f {[G(a - u) - G(-a))1(u,v) 

•'0 

+ [G(a) - G(-a + u)]j(-u, v)} du. (35) 

Using equation (35) and assuming uniform illumination, E a (x,y) = 1, 
two integrations are readily eliminated and equation (28) reduces to 



</) av = 4 exp (- fi 2 ) f [ (2a - u)(2b - v) 

Jo Jo 



■ exp [8 2 r(u, v)] cos p x u cos fi v v du dv. (36) 

By expanding the exponential function in a power series, equation 
(36) can bo divided into two parts corresponding to the coherent and 
incoherent parts of gain, as follows 

</> av - I. + I int (37) 

where 

/, = A'exp(-^)( s y a - asi y )' (38) 
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and 

h n c = 44 2 f f ~ —, ^ cos P:U cos W dv du - AI ine (39) 
where 

f«V(t*,»)l" 



M inc = 4 exp (-6') £ / / [2(6w + av) - uv] 

n - l "O "O 



n! 

• cos f} x u cos /3„z; rfw dv (40) 
and 4 = 4a b is the aperture area. 

The coherent part of the gain is the same as the antenna gain in 
the absence of surface deviations but multiplied by exp(— 8 2 ). This 
follows from equation (28) by expanding the exponential function 
which contains the correlation function in a power series. 

To obtain an estimate for the on-axis gain, equation (39) is eval- 
uated for fi x = p y = 0, and for a gaussian correlation function with 

r(u, v) = exp f - —~J^) ( 41 ) 

where c is the correlation distance. 
On-axis 

IUO, 0) = *Ac> £ ^- [l - ~-> f 1 + l) + -T-l • (42) 

In equation (42) terms of order exp[— n(2a/c) 2 ] and exp[— n(2b/c) a ] 
were neglected. 

By extending the limits of integrations in equation (39) to oo, the 
integration of the first part of this equation can be performed and 
gives equation (11). 

APPENDIX B 

On- Axis Gain for Circular Aperture Antennas 

For circular aperture antennas the on-axis gain for uniform il- 
lumination and a gaussian correlation function is obtained from 
equation (8) by expanding the exponential function and performing 
the integrations for the n = term, and the integration with respect 
to the azimuthal coordinates for the remaining terms, resulting in 

m , a > U .ffi mfi -f) + (*)'±i££ Im (43) 
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where D is aperture diameter, and 

I., - j°" (" csp [- Sfc^]^), dm dp , . (44) 

I = Modified Bessel function of order zero. 

The two integrations in equation (44) can be performed either with 
the aid of the Q(y, a„) function defined by 11 

Q(y, a„) = j^ exp (- ^^)l (xy)x dx (45) 

or by means of recently evaluated integrals of products of Bessel 
functions. 12 Let us use the former method. Let 

x = (2n) i p/c (46) 

y = (2n)Vc (47) 

a n = D/c(n/2)K (48) 

With equations (45) through (48), equation (44) can be written 

Ln = (jjg)' [" [1 - Q(y, a n )]y dy. (49) 

Integrating by parts results in 

/..-0f[l-« a ..a.), + /;^4 (50) 

The derivative in equation (50) can be expressed as 

|2 = a „ exp (_«li£) /lM . m 

Equation (51) is readily derived from (45) and the following integral 13 

f exp (-f/2)J.(xf)J,(yt)t clt = exp (- - + y )l.(xy) (52) 

where J, is a Bessel function of order v. Substituting equation (51) 
into (50) and integrating by parts yields the result 

Ln = ( C ^) 2 {f [1 - Q(fl* , a„) - exp (-aDlU)} 

+ 1 P exp [ ~* (a " + y2)] ij w< a »y )] dy ] • (53) 
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Using the relation 

— [yli(a n y)] = a n y I (a„y) (54) 

and the definition of the Q function, equation (45) , yields 

Ln = (^){o»[l " Q(an , On)] - f exp (-aDlMJ- (55) 

To evaluate Q(a„,a n ), the following relation, readily derived from 
(45) and (52) , is used. 



Q(a, 0) + QQ3, a) =2 + JT exp (-tf){f t [Jo(af)J {fif)]j *t. (56) 
Integrating (56) by parts and using (52), gives the known relation 

Q(y, x) + Q(x, y) = 1 + exp {-^^)h{xy). (57) 

With equations (55), (57), and (48), (44) is given by 



'-(TO—t-sOK-a+'-S? 



(58) 



The gain on axis (43) can therefore be written by using equation 
(58) as: 

0(0, 0) = (£) 2 [s 2 exp(-5 2 ) + (|) 2 5 2 exp(-5 2 ) ±Jfa [1 - Aj] 

(59) 



with 



A„ = 



-h(9l^fl+^fl]- « 



For n/2(D/c) 2 » 1, when the modified Bessel functions can be 
approximated by the first terms of the asymptotic series, A„ is then 
given by equation (15). 
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